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Abstract 

The most general situation of the reinforcement of a plate with multiple holes 
by several patches is considered. There is no restriction on the number and the 
location of the patches. Two types of the patch attachment are considered: only 
along the boundary of the patch or both along the boundary of the patch and the 
boundaries of the holes which this patch covers. The unattached boundaries of 
the holes may be loaded with given in-plane stresses. The mechanical problem 
is reduced to the system of singular integral equations which can be further 
reduced to the system of Fredholm equations. A new numerical procedure for the 
solution of the system of singular integral equations is proposed in this paper. It 
can be observed that this procedure provides a noticeable improvement from the 
one presented in Zemlyanova [§] , Zemlyanova and Silvestrov [9"] , To'] and allows 
to achieve a significantly better numerical convergence with less computational 
effort. 

Keywords: Reinforcement, patch repair, complex potentials, singular integral 
equations. 



1. Introduction 

Patch repair is one of the most common techniques used to strengthen thin 
constructions containing imperfections such as cracks and holes. There is a 
significant amount of scientific literature dedicated to the topic with different 
types of the plate-patch attachment being considered. In particular, the plate 
and the patches can be attached to each other using adhesive applied along parts 
of their surfaces (two-dimensional) or along some contours (one-dimensional), 
or only at certain point by using rivets. The detailed review of the literature 
on the bonded (two-dimensional) repair of the cracked plates can be seen in 
[3, 0, @ and many others. Reinforcement of plates with holes is studied less 
extensively. Engels et al. |3l have considered a bonded repair of an anisotropic 
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plate with a hole. Some problems on the reinforcement of a cracked plate with 
circular holes have been studied numerically in Q • The patch repair of a single 
hole in a thin plate reinforced by a single patch is considered for two different 
types of the plate-patch attachment in To the author's knowledge, a 

patch repair of multiple holes in a plate by multiple patches of arbitrary shapes 
has not been a subject of many studies with an exception of [loj . 

In this paper, a reinforcement of a plate with multiple holes with several 
patches is considered for two types of the plate-patch attachment: only along 
the boundary of the patch or along both the boundary of the patch and the 
boundaries of the holes covered by the patch. It is assumed that the patch is 
loaded with in-plane stresses applied to the free unattached boundaries of the 
holes and at infinity point of the plate. The location, the quantity and the shape 
of the patches can be arbitrary with the only condition that the boundaries of 
the patches and the holes are smooth and do not touch each other or intersect. 

The complex analysis methods are used to solve the problem. First, the 
Muskhelishvili's formulas [Hj are used to describe the stresses and the displace- 
ments in the plate and the patches through the complex potentials in the plate 
and the patches. The singular integral representations of the complex poten- 
tials are used then to reformulate the problem as a system of singular integral 
equations. These representations were first proposed by Savruk [7| to study the 
problems for cracked plates. The uniqueness of the solution of the resulting 
system of the singular integral equations can be proved in a manner similar to 

One of the main goals of this paper is to propose an alternative technique 
for a numerical solution of the derived system of singular integral equations. 
The method of mechanical quadratures applied for the solution of similar sys- 
tems in [§| , , [13] is difficult to implement for a large number of contours due 
to decreasing precision. The idea behind the method proposed here is in the 
approximation of the unknown functions by trigonometric polynomials. Sub- 
stituting these approximations into the system of singular integral equations 
reduces this system to the system of linear algebraic equations for the coeffi- 
cients of the trigonometric polynomials. Numerical examples show that even 
the systems of the relatively small order are sufficient to achieve a good accu- 
racy of the results. By considering numerical examples it can be shown that 
this method demonstrates a noticeably better convergence, larger independence 
from the symmetry of the construction and, hence, can be applied to a wider 
range of practical problems. 

2. Boundary conditions 

Consider an infinite thin elastic plate S weakened by several holes with the 
boundaries Lj, j = 1, . . . , n -I- r. Thin elastic patches 5*^, k = 1, . . . , n -f m, 
with the boundaries Tk are used to reinforce the plate S (fig. [T|). It is assumed 
that the patch Sk, k — 1, . . . ,n, covers the hole and is attached to the plate 
S both along the boundary of the hole Lk and the boundary of the patch F^. 
Other patches Sk, k — n + 1, . . . , n + m, arc attached to the plate only along 
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Figure 1: A plate S with holes and attached patches S^. 

their boundaries Tk and ean be located anywhere in the plate S, in particular, 
they may cover some of the holes Lj, j = n + 1, . . . ,n + r. Given in-plane 
stresses Pj{t) act on the free boundaries of the holes Lj, j = n + 1, . . . ,n + r. 
It is assumed that the lines Lj and are smooth Lyapunov curves and do not 
touch or intersect each other. 

The plate and the patches are homogeneous and isotropic. Their thicknesses, 
shear moduli and Poisson ratios are given by h, /i, v and hk, ^ik, v^. k = 
1, . . . , TO + n, correspondingly. The principal in-plane stresses af and are 
applied at infinity of the plate and act in the directions constituting the angles 
a and a + 7r/2 with the positive direction of the real axis. 

Assume that the patches and the plate are joined perfectly along the junction 
lines. This means that the displacements are equal both in the patch and the 
plate from both sides of the junction lines and the stress equilibrium conditions 
are satisfied on the junction lines. These conditions together with the given 
stresses on the free boundaries of the holes lead to the following set of boundary 
conditions: 

{u + iv)+{t) = {u + iv)l{t) = {u + iv)l{t), t€Lk, k = l,...,n, (2.1) 

h{an + iTn)~^{t) + hk{<Jn + iTn)t{t) = /lfc(<7„ -|- iTn)k{t), (2.2) 

t € Lk, k = l,...,n, 
{(Tn+iTn)t{t) =Pk{t), t € Lk, k = H + 1, . . . ,n + r, (2.3) 
{u + iv)+{t) = {u + iv)-{t) = {u + iv)+{t), t€Tk, k = l,...,n + m, (2.4) 

/l(cr„ -I- iT„) + (i) + hk{(Jn + iTn)t{t) = h{(Jn + iTn)~ it), (2.5) 

t e Tk, k = I, . . . ,n + m,, 
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where (m + iv){t) is the vector of the displacements at a point t of the plate 
S or the patch Sk, and ct„ and t„ are the tensile and the shear components of 
the stress vector acting on the tangent line to the curves Lfe or Tf^. Parameters 
without a subscript are related to the plate S; parameters with a subscript "fc" 
are related to the patch Sk- Here and henceforth, superscripts "+" and "— " 
denote the limit values of the stresses, the displacements and other parameters 
from the left-hand and the right-hand sides correspondingly of the lines Lk or 
Ffc. The direction of the curve Lk is chosen to be clockwise and of the curve Tk 
counterclockwise (fig. 1). 



3. Integral representations of complex potentials 

The stresses and the derivatives of the displacements in the plate S or in the 
patch Sk (k = 1, . . . , m + n) on any contour L' lying in S or in 5"^ correspond- 
ingly can be obtained from the Muskhelishvili's complex potentials ^(z), ^(2) 
by the following formulas [H]: 

K + iT„)(i) = $(t) +¥(!) + ^(t¥^ + *(!)), (3.6) 

at 

2^1(„ + ^v)(t) = K<S>{t) -m- f + ^ e 

For the patch Sk, all the parameters and the functions in these formulas should 
be taken with the subscript "fc" . In the case of plane stress take k = (3 — i^)/(l + 
v) and Kfc = (3 - i^fc)/(l -I- i^k), k ^ 1, . . . ,m + n. 

We will derive the integral representations of the complex potentials in the 
spirit of . The main idea is to divide the "plate-patches" system into separate 
components and introduce the unknown functions describing the jumps of the 
stresses or the derivatives of the displacements on each of the lines Lj and Tk- 
These unknown functions may be chosen in such a manner that some of the 
conditions (j2.1|) - ()2.5|) are satisfied automatically. 

First, consider the infinite plate S with the holes Lj which is subjected to 
the given in-plane stresses at infinity and on the lines Lj, j = n + 1, . . . , n + r. 
The stresses in the plate S have a jump discontinuity on the curve 7 = ^k=T^i^ 
lying in the interior of the plate, and the displacements are continuous on these 
curves: 

qit) = {{an + iT„)+(t) - {an + ir„)"(i)) /2, t G 7. (3.7) 

Extend formally the plate S to the full complex plane so that the stresses are 
continuous through the curve L = U"^[Lj and the derivatives of the displace- 
ments have a jump discontinuity along this line: 
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Then the complex potentials <J>(z) and ^'(z) can be taken in the form [7]: 

^- _^ ' 9'it)dt , (n + l)-^ [ q{t)dt 

k=n+l 



t - z 



7i-\-r 



z - 2^ Jl \ t-z {t-zy ) 



{k + I)-^ f / Kq{t)dt tq{t)dt 



z G S, 



r= ((7^ + cr^)/4, r' = (a^ -ai°°)e-2*"/2, 
2vr(l + k) 

where Zfc (fc = 7t, + 1, . . . , n + r) is the arbitrarily fixed point inside of the contour 
Lk- 

Formally extend each patch Sk, k — 1, . . . , n, to the full complex plane so 
that outside of the patch all the displacements and the stresses are equal to zero. 
The displacements in the patch Sk are continuous through the line Lk C Sk-, 
and the stresses have a jump discontinuity on this line; 

qkit) = ((a„ + iT„)+(t) - K + ZT„)-(i)) /2, t e Lk. (3.10) 

In the complex plane thus obtained from the patch 5*^ the jump of the vector of 
the stresses (cr„ + ir„)fe(t) acting on the tangent line to the contour is equal 
to (ct„ + iTn)^{t), and the jump of the vector of the displacements {u + iv){t) is 
equal to {u + iv)'^{t): 



g'kit) = T^^4(" + ^v)+{t), t e r,. (3.11) 

i[Kk + 1) dt 



According to the condition (|2.5p and the formula p.7p we have ((t„ +*T„)^(t) = 
—2d'j^^q{t), t G F/j, dk = h^h. Therefore, the complex potentials 5'fc(z) 
can be taken in the form 

m Ji^ t - z 2tt J-p \ Kk + 1 J t - z 



.T, / N (^fc + 1) ^ f / Kkqk{t)dt _ tqk{t)dt 



2 



1 

2^ 



2iKkq(t) \ dt ( 1 1 \ J, ^ \ tdi 



dkiKk + i)Jt~z V + 1 7(^-2;)^ 

z e Sk, fc = 1, . . . ,n. 
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Observe that the representations p.l2p can be also used if the patch Sk 
covers several holes and is attached both along its own boundary and all 
the boundaries of the holes which are covered by Sk- In this case the line Lk 
needs to be understood as the union of all the boundaries of the holes which are 
covered by Sk- 

The same representations (|3.12|) are also valid for the patches Sk, k — n + 
1, - . - ,n + m, assuming that qkit) = 0. 

Hence, the stressed state of the "plate-patches" system is described by the 
complex potentials (|3.9|) . p.l2[) which contain unknown functions g'{t), t £ 
Lj, j = l,.--,n + r; qj{t), t G i^, j = and q{t), g'kit), t G Tfc, 

k = 1, . . . , n + m. We will look for these functions in the class of functions 
satisfying the Holder condition on the corresponding curves Lj and Tk- This 
choice guarantees the existence of all the principal and the limit values of the 
integrals of the Cauchy type in the formulas p.9|) . p.l2p . 

4. The system of singular integral equations 

To find the unknown functions g' {t), qk{t), q{t), g'k(t) we have the conditions 
(|2.1l) - (|2.4p and (|3.1ip . Recall that the condition (j2.5p has been used to derive 
representations (|3.12l) and thus is satisfied automatically. Observe also that 
the condition appears from the way we extend each patch 5*^ to the full 

complex plane (that is, it guarantees that the displacements and the stresses 
outside of the patch 5*^ are equal to zero) . 

From the formulas (I3.6|) and the representations p.9p . (I3.12p . satisfying the 
conditions (|2.ip - (|2.4l) and (|3.1ip . we obtain the system of singular integral equa- 
tions on the closed curves Lj, j — I, - . - ,n + r, and Tk, k — 1, . - - ,n -\- m, with 
the unknown functions g'{t), qkit), q{t), g'k{t)'- 




j=n+l 



( 



t — Zj t — Zj dt \t — Zj (t 



nQj Qj dt / kQj 



tQj 




) 



) 
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(^fc + l)-l f ( 1 T-t dt\ ^ 

Iff 1 T-t di\ ,^ 

^Ai^)k Kv^t^y^tjt)''^^^''^ 

+ ^KlVT)l(^-(f^l)^' fc = l,---,n; (4.13) 

Iff! T-t di\ ,^ 
m J^\T-t f-tdt)^ ' 

(4.14) 

t G Lk, fc = 1, . . . ,n; 

TTdtJ^/^^ 27r Ji^\T-t f-td.tj''-' 

I [ ( I T-t dt\ , , 
+ 2;^iAF^--(f31)2^J^W'^- 

+ (-+^)'^ / f ^ - ,iT)dT 

ni J^\T-t f-tdt)^ ' 

=-r-.i.;f (..(^)-l(^.a;))-'<'). 

i G Lfe, = n + 1, . . . , n + r; (4-15) 
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- f - f '— - V f «Qj _ + — ( ^'^J ^ ^ 

dt _^—^_^\t — Zj i ^ Zj dt \t — Zj {i — ZjY ^ 



TT C?t JPfc 27r 7l VT — t T — t dt 

Iff 1 r-i di 



2tt Jl\ f-t (f - f)2 
(k + 1)^^ [ ( ^ K dt 



iri \ r — t T ^ t dt 

(k+ 1)^1 /■ 1 T-t lit 

f ~i ~ [f-tfH 



g'{T)dT 
q{T)dT 



q{T)dT 



2 TTl ^Lfc V'''^^ T — t dt 

{nk + iy^ f f ^ T-t di 



f ~t [f -tf dt 



qk{T)dT 



1 f / Kk I dt 



2tt \t — t T — t dt 
Iff 1 T -t dt 



9k{r)dr 



2tt Jy^^\ f - t (f - i)2 dt 
id'^^ f f Kk Kk dt 



9'k{r)dT 
q{T)dT 



TT{Kk + 1) JPfc V"^^^ T — t dt ^ 

t G Ffe, fc — 1, . . . , n + to; 

where qk{i) =0, fc = n + 1, . . . , n + m, and dk — hk/h, k — 1, . . . ,n + m. 

Observe, that the additional terms /l^. /rj. )'^^ 

have been included into the equations (I4.15P and (I4.16P correspondingly. Adding 
these terms does not change the conditions (|2.1l) - (|2.4p and (|3.11l) and provides 
that the following conditions are satisfied: 

/ 5'(r)dr = 0, fc = n + 1, . . . , n + r, (4.17) 

/ g(r)dr = 0, fc = l,...,n + TO. (4.18) 

From the physical viewpoint the condition (|4.17l) means that the displacements 
are single- valued in the plate along the contours Lk, k = n + 1, . . . ,n + r, while 
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the condition ([4.181) provides that the resultant force apphed to the patch Sk is 
equal to zero [8[ , [10| . 

It can be shown that the system (|4.13l) - (|4.16p has a unique solution. This is 
accomplished by reducing the system ()4.13|) - (|4.16p to the system of Fredholm 
equations and proving the uniqueness of the solution of the mechanical problem. 
This has been proved in details for the case of the plate with one hole reinforced 
by one patch in . The proof for the general case can be obtained using similar 
methods. 

5. Numerical procedure 

In this section we present an alternative technique for solving the system of 
singular integral equations (|4.13|) - (|4.16p . The method of mechanical quadratures 
presented to study the problems for cracked plates in Q and adopted for plates 
with holes reinforced by patches in [§[, [9|, [l3| provides a sufficient accuracy 
with a reasonable computational time for the systems of similar type in the 
cases when the number of contours Lj and Tk is small. However, the results 
become unsatisfactory both in the sense of accuracy and computational time if 
a larger number of contours is considered. Another observed disadvantage of 
the method of mechanical quadratures is that the results obtained using this 
method tend to deteriorate significantly with the loss of the symmetry in the 
construction. 

To overcome these disadvantages the following approach is proposed and 
shown to be computationally effective. First, parametrize each curve Lj and F^ 
with a parameter 9 G [0, 2tt]. Let tj = Tj{9) be a parametrization of the curves 
Lj, and = ffc(^) be a parametrization of the curves F^. Next, approximate 
the unknown functions g' (t), qk{t), t E Lj, and q{t), g'^it), t £ F^, by the partial 
sums of the complex Fourier series by the parameter 0: 

N 

m= — N 
N 

E ^r^™'' J = l,...,n; (5.19) 

m=-N 
N 

9W= E «o"e™'. ^er,. j = l,...,n + m; 

m=-N 
N 

E 5^e""^ tSF,, j = l,...,n + m, 

where g'^j , , , g™ are the unknown complex coefficients of the partial sums 
(|5.19p which need to be found. 
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Figure 2: A plate with two square holes reinforced by two square patches. 



Substituting the formulas (|5.19l) into the system (|4.13l) - (|4.16|) leads to the 
following equations: 



n+r / N N 




El E 9hiA^:kio)+ E 9^,Bj:,ie)\ 

N \ 



m=-N 
N 



(6)+ J2 <lTBTk+2n+,n+riO)] (5-20) 



m=-N 



N 



m=-N 



6'e[0,27r], j ^l,...,An + 2m + r, 

where N is the degree of the trigonometric polynomials (|5.19p . and the func- 
tions AJ^^.{e), B'^^iO) and fj{e) can be obtained from the system (|433l) - (|4J6l) . 
The integrals in these functions can be computed numerically by using any 
appropriate method. 

The equations (I5.20p . in general, can not be satisfied for all the values of 
the parameter 9 due to the approximation of the unknown functions by the 
truncated Fourier series (|5.19l) . Instead, we assume that the equations (|5.20p 
are satisfied only in the finite number of points Og = tt{2s — 1)/{2N + 1), 
s = 1, . . . , 2N + 1. This transforms (I5.20[) into the system of linear algebraic 
equations with real coefficients of the order 2(2iV -I- l)(4n + 2m -|- r). Thus, the 
approximate solution of the system of singular integral equations (|4.13p - (|4.16p 
can be obtained by solving the system of linear algebraic equations. 

While the precise investigation of the convergence of the solution of the 
resulting system of linear algebraic equations to the solution of the system of 
singular integral equations (|4.13p - (|4.16p is outside of the scope of this paper, 
the numerical results obtained by using this method and also comparison with 
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e 

Figure 3: The comparison of the approximations of the real and imaginary parts of the function 
9[iTi{9)) for AT = 10, TV = 15 and N = 20. 

the numerical results obtained by the method of mechanical quadratures in [loj 
suggests that even relatively small number of terms TV {N < 20 in most cases) 
is sufficient to obtain a satisfactory convergence, thus, reducing the problem to 
the solution of a relatively small system of linear algebraic equations. 

For instance, consider a plate with two holes Li and L2 in the shape of 
the squares with rounded corners given by the equations ti^2{0) = 0.45(e*^ + 
g-3ie^g-j -p 2 which is reinforced by two patches in the shape of the squares 
with rounded corners with the boundaries Fi and F2 given by the equations 
^1 2(6') = 0.7(e*'' + e-3'^/14) =F 1 (%• E]). Assume that the patches attached to 
the plate along both contours Li, Fi and L2, F2 correspondingly. The plate is 
loaded at infinity by the tensile stress acting in the direction constituting 
the angle a = tt/A with the positive direction of x-axis, and = 0. The 
mechanical parameters of the plate and of the patches are fj, = 60, 1/ = 0.4 and 
Ml = M2 = 40, vi ^ 1^2 ^ 0.3 correspondingly. The thicknesses of the patches 
and of the plate are equal h = hi — h2. The fig. [3] shows the approximations of 
the real (solid line) and imaginary (dashed line) parts of the unknown function 
gi{Ti{9)) plotted for different numbers N of terms of the trigonometric poly- 
nomials ([519]) : N = 10, N = 15 and TV = 20. It can be seen that a good 
approximation is obtained already for iV = 15. Similar results have been ob- 
served in other considered examples. The comparison of the results with the 
method of mechanical quadratures [lo| is presented on the fig. SI The computa- 
tion is made for the same mechanical and geometric parameters of the plate and 
the patches as in fig. [2l and crj^ — 1, af = 0, a = 0. The patches are joined 
with the plate only along their boundaries Fi and F2. The solid-diamond line 
corresponds to the results obtained by the method of trigonometric polynomials 
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Figure 4: The comparison of the approximation of the function Im g'^ (ri (9)) computed by the 
method of trigonometric polynomials and the method of mechanical quadratures [lOll . 

with = 20 presented here and the dash-dot Hne corresponds to the method of 
mechanical quadratures with 100 points on each of the hues [13] ■ It can be seen 
from the fig. [4] that a good comparison of the resuhs is observed. At the same 
time it needs to be stressed that the method of mechanical quadratures tends 
to lose its efficiency for a larger number of contours and that is why alternative 
numerical techniques should be used. 

The fig. [5] shows the graphs of the tensile and the shear stresses (t„ and t„ 
on the boundaries of the left hole Li and the left patch Fi for the construction 
shown on the fig. [5J The plate is loaded at infinity by the tensile stress a"^ = 1 
acting in the direction constituting the angle a — t:/A with the positive direction 
of the real axis, and = 0. The mechanical parameters of the plate and the 
patches are the same as in the previous example. Here and further, on the 
hole boundary ij, the graphs given by the solid line represent the stresses in 
the plate, the graphs given by the dash line represent the stresses in the patch 
from the outside of the hole boundary Lj , and the graphs given by the dash-dot 
line represent the stresses in the patch from the inside of the hole boundary Lj . 
Similarly, on the patch boundary Tj , the graphs given by the solid line represent 
the stresses in the plate from the inside of Tj , the graphs given by the dash line 
show the stresses in the plate from the outside of Tj, and finally, the graphs 
given by the dash-dot line represent stresses in the patch. 

The graphs of the dependence of the maximums of the stresses <t„ and t„ on 
the lines Li and Fi on the patch size characterized by the parameter r are shown 
on the fig. [6l The parameter r here is the radius of the circle circumscribed 
around each square patch. The proportions of the patch stay the same as in the 
previous example and only the size of the patch is changing. The graphs on the 
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Figure 6: The dependence of the stresses on the lines Li and Fi on the size of the patch. 
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Figure 8: The dependence of the stresses on the lines Fi and F2 on the angle a. 
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Figure 9: A plate with tliree circular holes reinforced by three square patches. 

lines L2 and r2 are the same due to the symmetry of the construction. It can 
be seen from the graphs that while the maximums of the stresses on the hole 
boundary Li decrease as the size of the patch increases, the maximums on the 
patch boundary Fi increase with r. 

Consider the same construction as on the fig. [2] where the left patch Si is 
attached to the plate S along both lines Li and Fi while the right patch S2 
is attached only along its boundary F2. All other mechanical and geometrical 
parameters are the same as in the previous example. The plate S is loaded 
at infinity with the tensile stress = 1 acting in the direction constituting 
the angle a with the positive direction of the real axis. The dependence of 
the maximums of the stresses cr„ and r„ on the line Li on the direction of the 
loading given by the angle a is shown on the fig. [T] Observe that the stresses 
(T„ and r„ are equal to zero on the line L2. The graphs of the stresses on the 
lines Fi and F2 are shown on the fig. [H It can be seen that the direction of 
loading makes a significant impact on the maximums of the stresses and, hence, 
on the possibility of the failure of the construction. Thus, each potential repair 
should be made with a consideration for the direction of the loading which is 
applied to the construction. 

The fig. [9] shows a plate with three circular holes each of the radius R — 0.5 
and reinforced with three square patches 5*1, S2, S3 with the rounded corners. 
The centers of the holes and the corresponding patches are located at the points 
zi = — 1 — i, Z2 — 1 — i and Z3 — —l + i. The boundaries of the patches are given 
by the equations ^j{e) = 0.675 ■ e*(/3-T/4)(-gie _^ e-^'^Q) + zj, j = 1,2,3. The 
patches are attached to the plate both along their boundaries and the boundaries 
of the holes they cover. The angle f3 here represents the orientation of the 
patches (fig. |9|). The plate is loaded at infinity with a normal stress a"^ = 1 
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acting in the direction of the real axis {a — 0). The mechanical parameters 
of the plate and of the patches are /i — 60, ly — 0.4 and = /J.2 = Ms = 40, 
j^]^ = j/2 = = 0.3, /i = /ii = /i2 = /13 correspondingly. 

The graphs of the dependence of the maximal stresses cr„ and t„ on the lines 
Li, L2 and L3 on the orientation of the patches /3 are shown on the fig. [TOl The 
corresponding graphs on the lines Fi, r2 and T3 are shown on the fig. [TT] It 
can be seen that the maximums of the stresses in the plate and patches depend 
strongly on the orientation of the patches with respect to the loading applied. 
The minimums of the stresses are achieved at /3 « 0°, and the maximal stresses 
achieved at /3 « 60°. There is a noticeable difference in the magnitude of the 
stresses acting in the construction depending on the orientation of the patches, 
which means that an improper repair can increase the stresses even compared to 
the unreinforced case and can potentially lead to a failure of the construction. 
Thus each potential repair should be carefully evaluated in order to reduce the 
maximal stresses in the construction. 

The numerical method proposed here for the solution of the systems of sin- 
gular integral equations of the type (|4.13|) - (I4.16|) can be easily extended to more 
general constructions. It also has been observed to provide more stable numer- 
ical results for a wider range of geometries of the holes and patches compared 
to the method of the mechanical quadratures [!],[§], , to be less dependent 
on the symmetry of the construction and to be more computationally efficient. 

6. Conclusions 

In this paper we extended a recently suggested technique for solving the 
problem of an infinite elastic plate containing several holes of arbitrary shapes 
reinforced by multiple patches of arbitrary shapes to the case when each patch 
can be attached to the plate by one of the two methods: either only along the 
boundary of the patch or both along the boundary of the patch and the bound- 
aries of the holes covered by the patch. The problem is reduced to the system 
of singular integral equations. A new and more efficient numerical approach 
is proposed for numerical solution of this system. The numerical examples for 
some particular cases of the plate and patch geometries are given. Obvious fu- 
ture developments of the method are the implementation of finite boundaries of 
the plate and inhomogeneously imperfect interfaces and also the optimization 
of the patch repair in accordance with a given plate geometry and an applied 
loading. 
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